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Abstract 

We study the variation of the hadron masses in the presence of external magnetic 
fields of strength of the order of the masses themselves. 

We identify the main factors affecting the lattice simulation results: 

• the boundary discontinuities for eB << 27r/i^a^. 

• the SU(6) choice of the hadron wave-function. 

We confirm qualitatively the earlier theoretical ansatz on the linear behaviour 
of the masses with the magnetic field and, as a by-product, we improve the lattice 
measurements of the nucleon magnetic moments. 

However our systematic and statistical errors preclude us from measuring the 
theoretically predicted field strength at which the proton becomes heavier than the 
neutron. 
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1 Introduction 



We study the variation of the masses of the proton and the neutron in the 
presence of a uniform external magnetic field using lattice QCD. 

The behavior of proton, neutron and electron in the presence of a constant 
magnetic field B was considered in [1]. 

For a sufficiently small magnetic field, it was assumed in [f ] that the depen- 
dence of the hadron masses on B is controlled by their anomalous magnetic 
moments which are known from experiment. The masses of proton and neu- 
tron can be easily computed using Landau's formula for a Dirac particle in a 
constant magnetic field: 

\eB\ eB , , 

E: = m„ - ^g^s (f .2) 

s G { — denotes the spin, m,p and m„ the masses in the absence of B. 
The formulae hold for zero momentum and the lowest Landau level. 

Note that the |ei?| coefficients (f/2mp and respectively 0) were considered 
non- anomalous . 

If one inserts the experimentally known anomalous magnetic moments, 
gn = —1.9 and g,p = 2.79, then for increasing magnetic field the lowest state 
of the neutron goes down steeper than the lowest state of the proton. Conse- 
quently since in nature the neutron is heavier than the proton there exists a 
value of the magnetic field Bq where the masses of the two particle are equal. 
Bq is a very large field of the order of lO^'^T 10~^GeV^ . For magnetic fields 
larger than Bp lO^^T the proton, neutron and positron acquire masses which 
allow the /3-decay of the proton. 

Magnetic fields of that order of magnitude may appear in certain astro- 
nomical phenomena [2, 3] and an unstable proton would be a very interesting 
phenomenon. 

However the magnetic field required is very large so that the assumption of 
pointlike particles does not necessarily hold and the effect of the strong forces 
has to be taken into account. 

Since there is no analytic nonperturbative approach in QCD it is worthwhile 
to treat the problem numerically by means of lattice QCD. 

Fxternal magnetic fields have already been used as a tool to compute g- 
factors of proton and neutron [4, 5]. In the present paper we study the spin 
up/down states of proton and neutron and compare them which the phe- 
nomenological approach (eq. (1.1) and (1.2)). In particular we compute the 
slopes (Fnergy/i?) for the magnetic fields of interest. This is harder then 
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the computation of the g-factors where it is sufficient to measure the mass 
difference of spin up and down states which is more stable against system- 
atic and statistical errors. 

We use a method similar to the one described in [4, 5]: The magnetic field 
in the 2;-direction is applied to a set of SU (3) configurations by multiplying 
(2.1) the links by appropriate U{1) phases (eq. (2.6) or 2.8)). In this way we 
define the gauge background for computing the propagator (2.10) of a quark 
with charge e^. These propagators are then used in order to compute the 
correlation functions for proton and neutron (2.12). 

The exponential time dependence of these correlators define the nucleon 
masses m in the magnetic field: 

1 li Git) , , 

"^ = -i^7^°gG(o) ^^-^^ 

We find a dependence of the results on the different possible choices of the 
magnetic field configuration ((2-6), (2.8)). We describe how this problem 
is related to the gauge invariant definition of the charged hadron correlation 
functions (2.14). In order to get a definite procedure we consider first the 'fiat' 
case where no QCD is present and tune the relevant parameters such that 
in the 'fiat' limit our procedure reproduces the continuum theoretical results 
correctly (Table I). In order to estimate the reliability of our results, and their 
sensitivity to the wave function choice we compare the SU (6) wave function 
for the neutron and proton correlators (2.12) with the {u^^Cu)d function [6] 
(Table II). While in the usual measurements there is no difference between the 
two definitions, in our context, the electric charges quantum numbers are as 
relevant as the color and spin ones in characterizing the wave function and 
significant differences are expected. This is discussed in detail in section 2. 
Section 3 presents the results: g-factors and slopes for proton and neutron 
(Table III, Figure 3). 

We work with two lattice sizes and find a dependency of the results on the 
volume V . A linear extrapolation in IjV (compare also [7]) leads to agreement 
with the phenomenological predictions (1.1) and (1.2). 

The precision which was achievable does not allow a definite answer to 
the question whether the proton becomes heavier than the neutron in strong 
magnetic fields. However this work maps out the elements which play a role 
in the lattice simulation of hadrons in magnetic fields and opens the way for 
more massive studies. 

In particular, our results are consistent with the main assumptions con- 
tained in the formulae (1.1) and (1.2): 

1. The presence of QCD field has a drastic effect on the ground state vibra- 
tion energy lowering the slope of (mp|-|-mpj) from (2/3-|-2/3-|-l/3) |i?|/mg 
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to \B\lmp and the slope of (m„| + m„j) from (2/3 + 1/3 + l/3)|i?|/mg 



2. The QCD field has only a small effect on the splitting (m^i — mpi)/mp 
leading to the small corrections to the anomalous magnetic moments of 
the independent quark model. 

3. Note that we have nothing to say about the QCD mass difference mp(0) — 
m„(0) in the absence of magnetic fields. This quantity is clearly totally 
beyond today's technology and we treat it as an external datum as the 
phenomenological treatment of [1] does. 

More detailed conclusions are given in section 4. 



In this section we describe the ingredients of the method and how to tune them 
properly. 

We work with a set of 30 quenched SU (3) configurations generated at 
/3 = 6.0 by a variant of the Cabbibo-Marinari algorithm [8]. The configura- 
tions are separated by 500 sweeps with 2000 sweeps thermalisation. Statistical 
independence has been checked by comparing the variance of various Wilson 
loops for different bins of configurations. We use lattice sizes of 16x16x8x8 
and 24 X 16 X 16 X 12 (t X x X J/ X 2;). 

We work with Wilson fermions at K = 0.1475 and 0.1525 for the small 
and K = 0.145 and 0.15 for the large lattice. The quark propagators have 
been computed by means of a preconditioned conjugate gradient algorithm 
[9]. Since this has to been done for various values for B and much more 
computer resources are needed compared to 'nonmagnetic' mass computations. 
We work with pointlike sources. 

The remainder of this section is organized as follows: 
Section 2.1 describes different realizations of the magnetic field and their impli- 
cations on the results. Section 2.2 deals with the choice for the wave function 
appearing in the hadron propagator. 

2.1 The magnetic field 

On the lattice, a uniform magnetic field B can be introduced by multiplying 
links Ufi{x) with an external abelian phase U^{x): 



to 0. 



2 The method 
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where x = (t, x, j/, z). 

In order to get a uniform magnetic field in z-direction one can set (ignoring 
for the moment the boundary conditions): 

U^{x) = exp(+ze7j/) 
Uy{x) = exp{-teS,x) 

Uf{i) = 1 ^2.2) 
Ufix) = 1 

In that way the fiux through a plaquette in the x-y plane is given by 

exp{ieBa^) = exp(z(7 + i^)) (2-3) 

which defines the magnetic field as 

a^eB = -f + C (2.4) 

On a finite lattice with periodic boundary conditions (at least in direction x 
and y) , in order to get an uniform magnetic field, B has to be a multiple of 
an integer n: 

o „ 27r 

a^eB = -j—j-n (2.5) 
To see that consider a possible realization for the magnetic field: 

U^{x) = exp{—ieBa^yLx) for x = — 1 

= 1 else 

U^{x) = exp{ieBa^{x - xo)) (2.6) 

U^{i) = 1 

Uf{x) = 1 

where x = . . . — I etc. 

Xq is parametrizing the Polyakov loop degree of freedom. 
According to (2.6) the plaquettes in the x-y plane are: 



P^{x) = exp{ieBa^{—LxLy + 1)) for x = — l,y = Ly — 1 
= exp[ieBa^) else 



(2.7) 



Thus the magnetic fiux is homogenous on all the x-y plaquettes with the 
possible exception of the corner — 1^ Ly — 1. In order that the corner 
plaquette will equal all the others, B has to be quantized according to (2.5). 

Unfortunately the lowest magnetic field fulfilling (2.5) is very large for 
reasonable sizes of the lattice: If we take for example = Ly = 8 and 
2GeV (/3 = 6.0) we get eB GeV^ (note that in practice the lowest 
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field corresponds to a down quark with electric charge = — |e). For a 
proton this corresponds to a mass shift 6m = fiB of the order of one GeV 
{fi = magnetic moment of the proton). 

Thus in order to apply magnetic fields of interest we can either 

- increase and Ly or 

- give up the quantization of i?, which means that the field is not homoge- 
nous anymore. 

Since increasing the lattice requires an amount of computer resources not 
available to us we have implemented two realizations of B which give up quan- 
tization. 

The first realization stems from the observation that only one plaquette 
of configuration (2.6) destroys the homogeneity if B is not quantized. If this 
'exceptional' plaquette is far from the source we may assume that its infiuence 
is not too strong. For that reason we place the source in the middle of the 
lattice Xs = L^J2, = Ly/2 (while the 'exceptional' plaquette resides at 
X = Lx — 1^ y = Ly — 1) and increase the lattice size in the x-y plane. We 
refer to this realization as '1-ep' type ('one exceptional plaquette'). 

The second realization has been used in reference [4]. It applies fixed 
boundary conditions in the x-direction: 

U^{t,x,y,z) = for X = Lx — 1 

= 1 else (2.8) 

Uy{t,x,y,z) = exp{ieBa^{x — Xq)) 

The plaquettes in the x-y-plane are given by 

P^{t,x,y,z) = ioTx = Lx-l 

= exp{ieBa^) else 

We refer to this realization as 'fixed' type. 

The freedom of choosing various Polyakov loop values was not exploited in 
[4]. However we will see that the choice of Xq is significant in order to estimate 
the systematic errors originating in the finite size of the lattice and especially 
in the large fields related with the breaking of the quantization condition. 

In order to estimate the systematic errors related to the lattice discretiza- 
tion, magnetic background representation and the wave function choices, we 
compare in the 'fiat' case (all SU{3) links set to one) our simulation results 
with predictions from the continuum. 

We define the quark propagator 

Gg{x,t;xs,ts) =< q{x,t)q{xs,ts) > (2.10) 
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(q=u,d for up and down quark) and the SU{6) wave function [10] (e.g. for a 
proton with spin up) 

+ 2ut(3)ut(l)Ji(2) 

+ 2ut(2)ut(3)JKl) 

- ut(l)ui(2)Jt(3) 

- ut(3)ui(l)Jt(2) (2.11) 

- ut(2)ui(3)Jt(l) 

- Ui(l)ut(2)Jt(3) 

- u,i{3)u^{l)d^{2) 

- Ui(2)ut(3)Jt(l) ) 

(where the antisymmetrisizing sum over the color degrees of freedom has been 
omitted). 

With these definitions the zero momentum propagator (e.g. for a proton 
with spin up) is given by 

Gp^it; Xs, ts) = J2< *pt(^. t)^pAxs, ts) > (2.12) 

X 

and we define the (time dependent) ratio of the g-f actors of proton and 
neutron by 

|- Gpt-Gpl 1 _ |- Gpt-Gp| 1 

Qn same but p ^ n 

In Table I we compare it with the independent quark model result which is 
—3/2 [11]. It shows the (time dependent) ratio for the two types of magnetic 
fields defined before and for two choices of Xq. The data show that: 

- for small fields [eBa^ 0.01) the continuum value of —3/2 is reproduced 
very well by all three types. 

- for larger fields the '1-ep' type especially for Xq = x^ is superior. 

- for large fields eBa!^ >sa 0.05 lattice discretization and nonlinear effects 
become important. 



As opposed to magnetic moment studies where only energy differences be- 
tween spin states are important, in the present work, the behavior of the single 
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magnetic field type 


eBa^ 


t 


'fixed' 


'1-ep' 


'1-ep' 






Xq = 


xo = 


Xq = Xg = 




3 


-1.47 


-1.47 


-1.50 




4 


-1.48 


-1.48 


-1.50 


0.015 


5 


-1.48 


-1.49 


-1.50 




6 


-1.49 


-1.50 


-1.51 




7 


-1.50 


-1.52 


-1.53 




3 


-1.24 


-1.24 


-1.47 




4 


-1.24 


-1.24 


-1.46 


0.06 


5 


-1.29 


-1.28 


-1.44 




6 


-1.36 


-1.35 


-1.43 




7 


-1.46 


-1.46 


-1.42 




3 


-1.03 


-1.22 


-1.42 




4 


-1.04 


-1.19 


-1.38 


0.12 


5 


-1.13 


-1.23 


-1.33 




6 


-1.20 


-1.28 


-1.30 




7 


-1.27 


-1.36 


-1.28 



Table I: no QCD, 16 x 16 x 8 x 8, K = 0.11: SU{6) (time dependent) ratio 
of the magnetic moments gp/gn (t) for different types of the magnetic field 
configuration. 

'proton' and 'neutron' states in an external magnetic field is of particular inter- 
est. The dependence of these states on B in the independent quark model can 
be computed by applying the Landau formula (1.1) separately to each quark 
contained in the SU (6) wave function. In particular we find that the proton 
spin up and the neutron spin down state should not depend in the "fiat" limit 
on B while the states with opposite spin should go up for increasing B. In 
the absence of QCD, the dominant effect is the presence of a contribution 
\eqB\/2mq from the zero-point vibration of each quark. 

If we compute these states for different choices of Xq we find a drastic 
difference between Xq = Xs = L^J2 and Xq = 0. 

Figures 1 and 2 show the effective masses nip j fit) = \og{G{t -\- 1)/G{t)) at 
t = 5 (which is representative for all time slices t) for the '1-ep' type and Xq = 
resp. Xq = Xs = Lx/2. While the choice Xq — X s reproduces the expected result 
the choice Xq = shows the protons state becoming heavier than the neutron 
states for sufficient large B. This may mislead to the conclusion that the 
conjectured effect - the proton getting heavier than the neutron - can be seen 
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- 






proton down -© — 
proton up -B— 

neutron down 

neutron up 






























0.02 0.04 0.06 0.08 0.1 0.12 

eBa"2 

Figure 1: no QCD, 16 X 16 X 16 X 8: effective masses for proton and neutron 
in '1-ep' field for Xq = 0. O = pj, □ = A = n^, X = n^-. see text. 

already in the absence of QCD. 

However this "effect" is just the result of a gauge dependent definition of 
the propagators with respect to the electromagnetic gauge transformations. It 
disappears when the appropriate gauge invariant propagators which take into 
account the electromagnetic phase-path corrections are used. 
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2.4 



2.3 



2.2 



2.1 



1.9 



1.8 



1.7 



1.6 







proton down -« — 






proton up -a— 






neutron down 






neutron up 

















0.02 



0.04 



0.06 
5Ba"2 



0.08 



0.1 



0.12 



''igure 2: no QCD, 16 X 16 X 16 X 8: effective masses for proton and neut 
ti '1-ep' field for Xq = x^. O = pj, □ = A = n^, X = n^. see text. 
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A hadron propagator invariant under U{1) gauge transformations can be 
constructed by including a product of U{1) links which lie along a path con- 
necting the source and sink. 

For zero momentum it reads 

G{t; Xs, ts) = J2< t(f , t)m{xs, ts) > (2.14) 

X 

where U denotes the product of U{1) links connecting x^t to Xg^tg: 

U := U{x,t; Xn.tn) . . .U{x,,t,; x,_i,t,_i) . . .U{xi,ti; x,,t,). (2.15) 

The exact choice of the path is somewhat arbitrary. However, this is immaterial 
for ideal QCD simulations where the mass of the lowest state extracted from 
the propagator (2.14) is unique if the wave function ^(x, t) has an overlap with 
the physical state of interest (and if one extracts the mass at sufficiently large 
t-t,). 

The experiments discussed so far have been done with U = 1. In order to 
take into account the phase we made experiments with the variants oi U 7^ 1. 
At the beginning we still discuss the case where QCD is switched off and we 
can compare with the continuum independent quark model predictions. As an 
example we choose for the path the shortest path between {xg^tg) and (x, t). 
The nontrivial contribution in (2.15) is then coming from 

Uy{x^t) = eyiY>[ieBa^[x — Xq)) (2.16) 

This corresponds to both types of the magnetic field '1-ep' (2.6) and 'fixed' 
(2.8) but with the x-links at the boundary omitted. In this way we get 

U = exp iBa^{x - Xo){y - y,) ^ (2-17) 

where the sum runs over all charges of the quarks contained in the hadron 
of interest. 

We find that this choice 

- does not change (almost) the behavior for 

- corrects the wrong behavior for = such that it becomes identical 
with the case Xq = x^ 

The results discussed so far imply that for the 'fiat' case the lattice repre- 
sents the continuum if we use the '1-ep' type and Xq = x^ or include the path 
in the propagator as described before. Note that due to the Polyakov loop 
(finite size) effects, the 2 choices could have given in principle different results. 

This choice has been used for all the results presented in the sequel. 
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2.2 Wave functions for proton and neutron 

In order to extract hadron masses from the time dependence of the correlation 
function the wave function for the particle of interest should have maximal 
overlap with the physical state. One way to achieve a better overlap is certainly 
to use non pointlike sources [12]. Due to the interaction of the magnetic field 
type with the position of the source we have considered the simplest case of 
pointlike sources. 

We have found that the overlap is much more sensitive to the choice of the 
wave function than it is for 'non-magnetic' mass simulations. 

In order to find a wave function which is suitable for the present work we 
were again guided by the case where QCD is switched off. 

It is well known that the independent quark model predicts a ratio for the 
g-factors of proton and neutron, gp/gn = —3/2 [11], which is remarkably close 
to the experimental value of —1.47. The value of —3/2 can be computed from 
the magnetic moments of the 'proton' and 'neutron' states given by the SU (6) 
wave function. 

In order to estimate the reliability of our numerical results and their sensi- 
tivity to the changes in the wave function we compared the SU (6) results with 
the ones obtained with the {u^^Cu)d wave function. 

In the absence of electric charge effects, the two wave functions are equiva- 
lent and {u^^Cu)d is the usually used wave function. In the presence of electric 
charge effects the SU{&) function is the correct one, while {u^^Cu)d gives the 
incorrect answer —5/4. 

We used these differences for checking the consistency and reliability of our 
numerical procedures. We have computed the (time dependent) ratio of the g- 
factors for proton and neutron (2.13) for the two choices of the wave function: 
[uC^f5u)d and SU{6). Table II shows the results obtained on a 16 X 16 X 8 X 8 
lattice from 20 quenched configurations at at /3 = 6.0, the '1-ep' type magnetic 
field with Xq = and K = 0.1475. For small magnetic fields eBa!^ < 0.02 the 
SU{6) data are close to the experimental value of —1.47. However the data 
for the [uC^f5u)d wave function are close to the —5/4 of the fiat case. They 
increase for large time-slices and thus indicate a mixing with higher states. 
This implies of course that the [uC^f5u)d wave function has a poorer overlap 
with the physical state than the SU (6) wave function. From the table we also 
can see that for larger fields eBa^ > (9(0.01) the ratio decreases with increasing 
magnetic field. This has already been seen in the fiat case (Table I). 
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9p/9n 




eBa^ 


t 




SU{6) 




3 


-1.22 


-1.48 




4 


— 1.23 


— 1.52 


0.01 


5 


-1.26 


-1.52 




6 


-1.26 


-1.50 




7 


-1.29 


-1.54 




3 


-1.21 


-1.48 




4 


— 1.23 


— 1.52 


0.02 


5 


-1.27 


-1.52 




6 


-1.26 


-1.50 




7 


-1.27 


-1.53 




3 


-1.18 


-1.45 




4 


-1.18 


-1.48 


0.08 


5 


-1.24 


-1.48 




6 


-1.18 


-1.45 




7 


-1.07 


-1.38 




3 


-1.14 


-1.41 




4 


-1.10 


-1.38 


0.15 


5 


-1.10 


-1.33 




6 


-0.99 


-1.28 




7 


-0.79 


-1.10 



Table II: QCD, 16 x 16 x 8 x 8, K = 0.1475, '1-ep', xq = x,: (time dependent) 
ratio of the magnetic moments gp/gn (t) for the {uC^su)d and the SU{&) wave 
function. 



12 



3 Results 



In this section we present the data obtained from the numerical simulations. 
In particular we have computed the slopes 

slope(B) = (m(B) - m(0)) ^^^^^ , (3.1) 

eB 

the g-factors for proton and neutron and the ratio of their magnetic moments. 

Table /// shows the numerical results in comparison with the experimental 
data. The upper two parts of the table contain the data for an order of magni- 
tude for the magnetic field, where - according to the mechanism described in 
the introduction - the proton should become heavier than the neutron. The 
lower third of the table contains the data for a larger magnetic field, eB = 0.06 
resp. eB = 0.08 (These data stem from 20 configurations). The results for 
the different values of B are the same within the error bars. In the previous 
section we have already seen that for larger magnetic fields our lattice does 
not reproduce the precise continuum results even in the absence of QCD. 

The data stem from a 16 X 16 X 8 X 8 lattice at K = 0.1475 and 0.1525 
and from a 24 X 16 X 16 X 12 lattice at K = 0.145 and 0.15. The data differ 
for these two lattice sizes. In particular the spin up slopes differ by amounts 
which can hardly be explained by the error bars. If we extrapolate the data 
to infinite volume by assuming a a = ao + 1/Vai behavior we get results 
which are remarkably close to the experimental data (Table III columns 2 
and 3). Moreover this extrapolation brings agreement with the proton-decay 
scenario of [1]. To see that, compare in Table III, column 3 the infinite volume 
extrapolation of energy/B slope of the lowest nucleon states: while for the 
neutron the energy decreases with slope —2.1 ( neutron J, state in line 1) for 
the proton energy decreases only with slope —1.9 (proton | state in line 5). Still 
we refrain from considering these results as definite and reliable confirmation 
of the theoretical prediction of the proton decay [1]. For such a confirmation 
larger lattices with better statistics will be necessary. 

The statistical errors have been computed by dividing the configurations 
into bins and comparing the results for each of them. The systematic errors 
are hard to estimate. The most significant one certainly stems from the use of 
pointlike sources. They make it difficult to avoid the systematic overestimation 
of the masses caused by higher excitations of the spectrum. We have estimated 
this effect by extracting the masses from the correlation function for different 
windows of time-slices. 

The ratio of the g-factors does not suffer from the former effects, since they 
consist of ratios of masses. As already presented in Table II the ratios (2.13) 
are very stable for a quite large numbers of time slices. 
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phenom./ 
exper. 



V 



oo 



K = 0.145 
eB = 0.03 

24 X 16 X 16 X 12 



K = 0.1475 
eB=0.02 

16 X 16 X 8 X 8 



neutron J, slope 
neutron | slope 
neutron g-factor 



-1.9 
+ 1.9 
-1.9 



-2.1 
+ 1.3 
-1.7 



-2.1(4) 
+ 1.9(3) 
-2.0(4) 



-2.2(4) 
+3.2(2) 
-2.6(4) 



proton I slope 
proton I slope 
proton g-factor gp 



+3.8 
-1.8 
+2.8 



+3.5 
-1.9 
+2.7 



+3.4(4) 
-2.9(4) 
+3.1(4) 



+3.3(4) 
-4.8(2) 
+4.1(5) 



9p/9n 



-1.47 



-1.50 



-1.51(3) 



-1.53(4) 



phenom./ 
exper. 



V 



oo 



K = 0.15 
eB = 0.03 

24 X 16 X 16 X 12 



K = 0.1525 
eB=0.02 

16 X 16 X 8 X 8 



neutron J, slope 
neutron | slope 
neutron g-factor 



-1.9 
+ 1.9 
-1.9 



-1.8 
+ 1.5 
-1.7 



-2.4(4) 
+2.0(4) 
-2.2(4) 



-3.5(8 
+3.0(7 
-3.3(8 



proton I slope 
proton I slope 
proton g-factor gp 



+3.8 
-1.8 
+2.8 



+3.5 
-1.1 
+2.3 



+3.8(4) 
-3.2(4) 
+3.5(4) 



+4.5(8 
-7.4(8 
+6.0(8 



9p/9n 



-1.47 



-1.44 



-1.53(8) 



-1.7(2 



phenom./ 
exper. 



V 



oo 



K = 0.145 
eB = 0.06 

24 X 16 X 16 X 12 



K = 0.1475 
eB=0.08 

16 X 16 X 8 X 8 



neutron J, slope 
neutron | slope 
neutron g-factor 



-1.9 
+ 1.9 
-1.9 



-2.2 
+ 1.3 
-1.7 



-2.5(5) 
+ 1.5(3) 
-2.0(5) 



-3.1(6) 
+2.0(3) 
-2.6(6) 



proton I slope 
proton I slope 
proton g-factor gp 

9p/9n 



3.8 
-1.8 

2.8 
-1.47 



+3.3 
-1.7 
+2.5 
-1.48 



+3.2(4) 
-2.8(4) 
+3.0(4) 
-1.50(4) 



+3.0(6) 
-4.9(3) 
+4.0(6) 
-1.54(5) 



Table III: slopes of proton and neutron states, g-factors and their ratios. 
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Figure 3 shows the proton and hadron masses as a function of the magnetic 
field for the 24 X 16 X 16 X 12 lattice and K = 0.145. 



1.55 



1.45 



1.35 



proton down 
proton up -B- 
neutron down x 
neutron up 




0.08 



Figure 3: QCD, 24 X 16 X 16 X 12: masses for proton and neutron for K = 0.145. 
For better readabihty error bars are shown only for one state. 
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4 Conclusions and Outlook 



In the present work we showed which ingredients are important in order to 
study the variation of hadron masses in an external magnetic field. 

In particular we pointed out the factors infiuencing the lattice results: 

- the choice of the (SU(6)) wave function. 

- the location and the boundary discontinuity of the magnetic field. 

- the path-phase correction which makes the charged-particles propagators 
QED-gauge-invariant. 

- the location of the source 

We computed the slopes (mass/magnetic field) for the lowest spin up and 
down states of proton and neutron and compared them with the phenomeno- 
logical approach of [1]. 

We worked at very large magnetic fields which are expected to drive the 
proton to a larger mass than the neutron. We found the remarkable result that 
-within the available numerical precision- the assumptions leading to (1.1) and 
(1.2) hold for these very large fields. However in order to decide whether the 
proton becomes heavier than the neutron more precision is necessary. 

Further investigations thus should include: 

- a carefull study of the finite size effects 

- better statistics 

- smeared sources adapted to the charged particles problem 

- extrapolation to the chiral limit. 
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